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Abstract 

We investigate the critical points of the near field intensity in the vicinity of dielectric bodies 
illuminated by the incident light. It is shown that in the electrostatic approximation there are no 
local maxima of the intensity outside dielectric surfaces and the only possible critical points are 
either local minima or saddle points. Using the boundary charge method we investigate numerically 
the field distribution around star-like sets of prolate spheroids. The field enhancement is shown to 
achieve a value of several hundreds at critical points outside the surfaces of spheroids. 
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I. INTRODUCTION 



Progress in near-field optics is caused by development of both experimental technique 
and mathematical methods of calculations of the near-field structures. Near-field microscopy 
has been successfully used to overcome the diffraction limit and to achieve subwavelength 
resolution imaging of various surface structures |l( . Of particular interest is a phenomenon 
of localization and enhancement of the EM near field by metallic nanoparticles and clusters 
Q]. For the explanation of the phenomena of the Surface Enhanced Raman Scattering 
of particular interest are large local fields induced in close proximity or at the surfaces of 
nanoparticles [s,^. Mechanism of extraordinary Raman enhancement is still not completely 
understood. However, according to electromagnetic mechanism fields localized on surfaces 
of nanoparticles are responsible for extraordinary amplification of the Raman scattering. 
According to electromagnetic mechanism P, 0, Q] the cross-section of Raman scattering is 
proportional to the fourth power of the local field. Therefore, even moderate enhancement 
of the local electromagnetic field results in the extraordinary amplification of the Raman 
scattering. Enhancement of the Raman cross section up to ten power fourteen was reported 



in literature thus allowing for single molecule Raman spectroscopy and detection 

Exact solutions to Maxwell's equations are known only for special geometries such as 
spheres, spheroids, or infinite cylinder, so approximate methods are in general required. 
Comprehensive review of methods of anal ysis used in near field optics theory may be found 
in review papers and in monographs ll^l 14]. Expansion of solution on multipolar 



eigenfunctions has been applied to many near-field optical problems in early studies jl^ . 
Several approximate methods have been used for calculation of spatial distribution of near- 
field induced by external electromagnetic wave around nanoparticles of arbitrary geometrical 
shape in electrostatic approximation, which is applicable when size of nanoparticles is much 
less than wavelength of EM field. Among them most popular are discrete-dipole approxima- 
tion (DDA) method or coupled dipole approximation (CDA) method Weak point 
of DDA or CDA methods is that they can not be applied for particles with large (exceeding 
2) refractive index because of its low accuracy and can not describe accurately fine structure 
with large field gradients. 

Another approximate method commonly used for calculating distribution and enhance- 

n 

ment of local fields is finite difference time domain (FDTD) method 18]. This method 
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was applied to calculate local field distribution in a junction between two silver nanoparti- 
cles [19]. This method also has disadvantage that it can not describe accurately large field 
gradients. 

For our purposes we found out that most suitable is Boundary Chage Method (BCM 



when Maxwell's equations are solved numerically using integral equation formulation 
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221 ] . Below we will show that by using this method one can more accurately calculate fine 



structure of local field distribution in space between plasmon nanoparticles and on surfaces. 

When metallic or dielectric body is exposed to the external electric field, the maximum 
field enhancement is achieved at the points of its surface with maximal curvature. Further 
field amplification is obtained by tuning the frequency of the external field towards the 
frequency of an appropriate surface plasmon oscillation. This is illustrated by Fig. ^ where 
the dependence of the frequency of the dipole plasmon resonance (a) and the maximal field 
enhancement (b) on the aspect ratio a/b (a > b) of a single prolate spheroid is depicted using 
the well-known solution of the Poisson equation in elliptic coordinates [^J . Here the external 
field is parallel to the large axis of the spheroid. In plotting Fig. a simple approximation 
for the dielectric function of silver, e{oj) = 6.4 — 100/c<j(u; + 0.06i) , where the frequency, u, 
is measured in eV, was used Q] • As is readily seen from Fig. ^ the field enhancement for 
a sufficiently long body may achieve several thousands. Also, increasing aspect ratio of a 
spheroidal particle results in a significant red shift of the plasmon resonance. 

Combining particles of various size and shape one may attempt to construct a nanolens 
with the region of the maximum field intensity (a hot spot or a nanofocus) located in a gap 
between nanoparticles apart from their surfaces. Recent computations implementing the 
multiple spectral expansion method j^] demonstrated the emergence of such a nanofocus 
in a system of self-similar chain of metallic spheres of different size. Such a possibility is 
very attractive since it allows for remote detection without the mechanical contact between 
a molecule and a nanoparticle. 

The main purpose of the present paper is to investigate the field distribution near hot 
spots in sets of plasmon nanoparticles. In section |HJ a simple classification of critical points 
of the field intensity is given. It is shown that in electrostatic approximation when one can 
neglect retardation effects absolute maxima of the field intensity are impossible in space out- 
side surfaces of nanoparticles and the only possible critical points outside surfaces are either 
absolute minima or saddle ones. The boundary charge method is explained in section IHII 
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FIG. 1: The plasmon frequency (a) and the maximal field enhancement (b) versus the aspect ratio 
of a prolate spheroid. 

and the results of numeric calculations of the near field distribution produced by symmetric 
star-like sets of prolate spheroids are presented in section ITVl Finally, in section IS the fine 
structure of the field distribution around a self-similar chain composed of metallic spheres 
[20! ] is investigated and it is shown that maximum field intensity is achieved on surfaces of 
nanoparticles. A brief summary is given in concluding section IVI1 

II. FIELD DISTRIBUTION NEAR CRITICAL POINTS 

When dimensions of nanoparticles are much less than the wavelength of an incident 
electromagnetic field, one can neglect retardation effects and use the electrostatic approx- 
imation to calculate the field distribution. Thus in the near-field region, the electric field, 
E(r) = — Vy?(r), satisfies Poisson's equation, 

Ve(u},r)Vp(r) = 0, (1) 
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with the boundary condition 

E (r -> oo) -> £ in , (2) 

where e(u;, r) is the complex dielectric permittivity and i£ m is the electric field amplitude of 
the incident electromagnetic wave. In traditional electrostatic problems, when the dielectric 
permittivity is independent of u, both the external field and the potential are real-valued 
functions. In application to the near-field optics, both the electric field of the incident 
wave, E m , and the potential, <p(r), may be complex- valued due to the possible elliptical 
polarization of the incident wave and due to the energy losses inside nanoparticles. 

Prior to discussing the results of numeric solution to Eq. (JTJ), we give the qualitative 
description of hot spots. To be more rigorous, we define a hot spot as a point where the 
field intensity distribution has a maximum in certain directions. In other words, we are 
interested in the critical points with zero gradient, V/(r) = 0, where the field intensity is 
I(r) = |-E(r)| 2 = |Vy?(r)| 2 . The question is what kind of critical points are possible apart 
from the surface. 

To answer this question suppose that the critical point is at the origin, r = 0. In a 
vicinity of this point, the solution to Poisson equation (Q) may be written as a sum 

tp(r) = y?o - E ■ r + <p 2 (r) + ip 3 (r) + . . . , (3) 

where Eq is the electric field at the origin and (p n {f) are harmonic polynomials of the order 
n. In particular, the quadratic and cubic parts of the expansion are 



ip 2 (r) = A 2 {3(ni • r)(n 2 ■ r) - (m • n 2 )r 2 } , (4) 
y? 3 (r) = A 3 {5(mi • r)(m 2 • r)(m 3 ■ r) (5) 
- r 2 [(mi ■ r)(m 2 • ra 3 ) + (m 2 • r)(m 1 ■ ra 3 ) + (ra 3 ■ r)(rrii ■ m 2 )}} 

where A 2t3 are arbitrary constants and ni i2 and mi i2j3 are arbitrary unit vectors. The field 
intensity near the origin looks like 

J(r) = El - 2Re (E* ■ V<^ 2 (r)) + | V V2 {r)\ 2 - 2Re {E* ■ V<^ 3 (r)) + . . . (6) 

Let us consider first the case of linear polarized incident wave; then E in in Eq. (J2J) is a 
real vector. If the losses inside dielectric particles are negligible, Ime(w) = 0, the resulting 
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solution to the Poisson equation is a real-valued function. According to Eq. © the condition 
of zero intensity gradient at r = is Eq ■ Vy?2 = 0. It is a matter of simple algebra to check 
that this yields either to E = or to A 2 = in Eq. (J1J). The first possibility, E = 0, 
corresponds to the absolute minimum of the field intensity, 1(0) = 0; this field configuration 
is used in the well-known Paul traps. The second possibility means that f2( r ) — in 
Eqs. (0BJ), hence the quadratic part of the expansion in Eq. © is —2E ■ Vy^r). The latter 
expression is a harmonic function. According to the well-known properties of harmonic 



functions [2^ it cannot take maximal or minimal values inside its domain. The only possible 
critical points are the saddle ones. Thus, with a linear polarization of an incident wave and 
negligible losses the field intensity apart from the surface of dielectric bodies may be either 
zero at certain points or it may exhibit a saddle point. 

In the case of complex-valued fields the analysis is similar but more tedious. The zero- 
gradient condition, V/(r) = 0, does not necessarily means that <f>2( r ) in Eq. Q vanish. In 
order to classify the critical points of the field intensity we have to investigate the eigenvalues 
of the quadratic form J 2 = | V</?2 (^) | 2 — 2Re (Eq ■ Vips(r)) appearing in Eq. (jHJ). Since all 
coefficients and vectors entering in Eqs. ()4J5j ) are generally complex, the quadratic form I2 
depends on 12 complex numbers. We have studied the signature of I 2 using Mathematica 
5.0 symbolic software. It was found that three eigenvalues are never negative all at once, 
that is, the field intensity © never takes maximal value apart from dielectric surfaces. The 
signature of the quadratic form, J 2 , may be either (+,—,—), (+,+,—) or (+,+,+). The 
latter possibility corresponds to the local minimum of the field intensity that may now exist 
even if Eq 7^ 0. Although such a "cold spot" is hardly of use in the near-field optics, it is 
worth writing down a corresponding example of the field distribution: 

<p(r) = -x + 2y + z + 2ixy + ^x 2 y - ^y 3 + ^x 2 z + ^xyz - ^y 2 z + ^yz 2 (7) 

Thus the most important conclusion from this analysis is that the field intensity cannot 
take maximal values apart from the dielectric surfaces. The spatial distribution of the 
intensity may exhibit either saddle points or minima. Below we demonstrate several numeric 
examples with both types of saddle points. 
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III. BOUNDARY CHARGE METHOD 



For the case of materially equal dielectric particles the problem is conveniently reduced 
;o a set of integral equations for the surface charge density, a, induced at the nanoparticles 



23|. Compared 



2l| or, equivalently, for the surface distribution of the normal electric field 
to various direct methods of the numeric solution of Eq. (0), the integral equation approach 
allows to reduce dimensionality of a problem and to achieve reasonable accuracy at smaller 
computer cost. Here, we use the boundary charge method (BCM) based on the integral 



equations in the form derived in 



2l|: 



where 



and 



A(u)a(s) = -n(s) ■ E in (s) + J ds'F(s, s')a(s'), (8) 



. . 1 + e{lS) , , 



n(s) ■ (s — s' 



F(s,s') = V v /|3 ' ■ (10) 

\s — s | d 

The integration in Eq. (JHJ) is carried over the multiply connected surface of all particles, 
vectors s and s' belong to the surface, and n(s) stands for the surface unit normal directed 
towards vacuum at s. Once Eq. (jHJ) is solved and the surface charge distribution, a(s), is 
found one can calculate the induced electric field at the arbitrary point outside the surface 
using a simple formula 

E(r) = E in - V r / ds t^-t. (11) 
J \r-s\ 

Eq. (jHJ) was approximated by a set of linear equations with the help of an appropriate 
triangulation of the surfaces. Our computer facilities allowed us to use up to 5000 triangu- 
lation points. As a test problem, we investigated the dipolar plasmon resonance for a single 
prolate spheroid. Comparing the numeric solution with the well-known analytic expressions 
(see fig. it was found that with the maximum affordable number of triangulation points 
the accuracy within 10% is achieved for the spheroid aspect ratio not exceeding 6:1. 
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FIG. 2: Sets of prolate spheroids used in computations. 

IV. RESULTS OF COMPUTATIONS FOR STAR-LIKE SETS OF NANOPARTI- 
CLES 

Now we turn to the discussion of the numeric solutions to integral equation (jHJ). We have 
investigated the field structure near the sets of prolate spheroids (Fig. |2J). For the examples 
discussed below, the aspect ratio of all spheroids is 4 : 1, their centers are at r = 4.5 from 
the origin. The complex dielectric permittivity, e(u), corresponds to silver. 

Fig. |21 shows the dependence of the field enhancement, / = |-E7|/|I2 m |, at the origin on the 
frequency of the incident wave for the configuration depicted in Fig. |2k- The incident wave 
propagating along the z axis is linear polarized with the electric field vector along the y axis, 
E vn = (0, 1,0). The field distribution along y axis is also shown. The electric field here is 
maximal at the surfaces of spheroids. At the origin, there is the maximum in zx plane and 
the minimum in y direction. The signature of the saddle point at r = is (— , +, — ). Notice 
that the maximum field enhancement for a single 4 : 1 prolate spheroid is about / ~ 300 
(fig. [TJ, that is, the maximum field at the origin is approximately half as much. With the 
circular polarization of the incident wave, there are just minor changes in the field structure. 

More complicated behavior is observed for the star composed of three spheroids (Fig. Et») • 
Fig. shows the field versus the frequency for the circular polarization of the incident 
wave, E in = (1, i, 0) /V2. The plasmon resonance frequency (2.42 eV) in this case is a little 
red-shifted compared to two spheroids (2.44 eV). The field structure, however, is entirely 
different. The field distribution along the y axis is depicted in Fig.HJ). The two-dimensional 
contour plot of the field intensity is shown in Fig. where the the electric field vectors are 
also plotted. 
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FIG. 3: (a) — The field at the origin versus frequency of the incident wave, (b) the field distribution 
along the y axis (oj = 2.44 eV) for the set shown in Fig. 



Instead of a single saddle point between two spheroids now there are four critical points. 
The minimum in the xy plane (marked with the x sign in Fig. |5J) is situated at the origin. 
Three others (O ) are maximums in the radial direction and minimums in the azimuthal 
direction. The field enhancement at these points is about \E\ ~ 150. Fig. 03 also shows the 
real and imaginary parts of the electric field vector. As is readily seen, the ellipticity of the 
wave is essentially nonuniform, i.e., this set of spheroids acts as a polarizer. 

In contrast with two spheroids, the field structure in the 3-star is sensitive to the polar- 
ization of the incident wave. Fig. |U] shows the intensity distribution and the electric field 
for the linear polarization of the incident wave, E m = (0,1,0). Now there are only two 
(+, — , — ) saddle points with the field enhancement \E\ ~ 130. Besides, the (+, +, — ) point 
moves to the lower part of the figure. 

Similar behavior was observed for the structures consisting of larger number of prolate 
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FIG. 4: a — Field amplitude versus frequency for the circular polarization of the incident wave and 
configuration depicted in Fig. Et". Solid line — r = (0, —0.45, 0), dashed line — r = (0, —0.45, 0). 
b — Field distribution along y axis. 



spheroids. Fig. shows the plasmon resonance for the circular polarization of the incident 
wave. The field distribution along a line coming through the gap between two adjacent 
spheroids is depicted in fig. [7jp. The two-dimensional plot (Fig. |HJ) demonstrates a single 
(+, +, — ) saddle point and four (+, — , — ) saddle points. 



V. FIELD ENHANCEMENT BY A SELF-SIMILAR CHAIN OF SPHERES 

As we have already mentioned, it was recently proposed to achieve larger field enhance- 
ment implementing self-similar sets of nanoparticles |2fjj. The idea is based on the fact that 
in the near-field region the field enhancement is independent of absolute size of particles. 
Therefore, if there two spheres of different size separated by a small gap, then the resulting 
field enhancement at smaller sphere should be about f 2 , where / is the enhancement factor 
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FIG. 5: Field distribution in xy plane around three spheroids (Fig. |2p) with circular polarization 
of the incident wave, x corresponds to (+, +, — ) saddle point, Q corresponds to (+, — , — ) saddle 
points. Green arrows — HeE, black arrows — ImE. 
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FIG. 6: Same as in Fig. |S] but with linear polarization of the incident wave. 

for a single isolated sphere, which depends on its material only. Combining m spheres of 
reducing radii one may expect to achieve the electric field at the surface of the smallest 
sphere enhanced by a factor of f m . 

We used BCM to investigate a combination of three silver spheres studied in jgoj]. The 
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FIG. 7: a — Field amplitude versus frequency for the circular polarization of the incident wave for 
the set shown in fig. b — Field distribution along the y = x line, c — Field distribution along 
the x axis. 

radii of spheres are R = 1, 1/3, 1/9, their centers are at z = 0, 1.53, 2.0. The electric 
field of the incident wave is parallel to the z axis. The field distribution along the z axis is 
depicted in Fig. El As is readily seen from the figure, the field at the surface of the smallest 
sphere is really enhanced by a factor of 600. However, the maximal field is achieved rather 
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FIG. 8: Same as in Fig. El for the 4-star set (Fig. Et). 

at the surface of the sphere then in a gap between spheres (see Fig. ^jp, where the magnified 
part of Fig. is shown). 



VI. CONCLUSION 



It is shown that in the nonretarded approximation the absolute maximum of the field in- 
tensity distribution around dielectric bodies illuminated by an external source is achieved at 
the surface. The only possible critical points of the intensity distribution outside surfaces are 
either saddle or minimal points. The computed field distribution around the sets of prolate 
spheroids having the aspect ratio factor up to 4 exhibits the maximum field enhancement on 
the order of 300 at the points of the surface with maximal curvature. It is worthy to note, 
that the field enhancement provided by a single silver sphere is about 30 However, in 
star like set of nanoparticles the field enhancement at the saddle points ("hot spots") outside 
the surfaces is also sufficiently large, / ~ 150. even for the aspect ratio 4. Since the Raman 
scattering is proportional to the fourth power of the electric field, the corresponding Raman 
cross-section is enhanced by a factor of 5 10 s . Preliminary computations with the increased 
aspect ratio and the fixed distance between spheroids demonstrated, first, the significant 
red shift of the plasmon resonance. Second, the field enhancement at the saddle points is 
approximately half the enhancement at the tip of an single spheroid with the same aspect 
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FIG. 9: Field distribution along a chain composed of three silver spheres, to = 3.35eV. Heavy lines 
mark the position of the spheres. 



ratio, that means, that field amplification at central points may achieve a value of several 
thousands for the aspect ratio of spheroids in a range 20-25 according to scaling law for 
isolated spheroidal shape particle presented on Fig. 1. That means that one could expect 
amplification of Raman cross section for molecule placed in the focus of star-like nanolens 
to achieve as much as 10 14 or even more. This may open possibility of the single molecule 
Raman spectroscopy in the scanning mode in an AFM type device with nanolense placed 
on the tip of AFM. 

The polarization state of the incident electromagnetic field is also of importance. Quali- 
tatively, the circular polarized wave may be represented as a sum of linear polarized waves 
with different directions of the electric field and appropriate phase shifts. In the symmetric 
configurations like those shown in Fig. |2l each wave excites its own spheroid. The field 
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distributions depicted in Figs. 14171 appear as a result of interference. 
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